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EXECUTIVE  SUMMARY 

This  report  details  a  feasibility  study  investigating  the  performance  of  a  nano-communication  system. 
The  study  is  based  on  analyzing  the  suitability  of  carbon  nanotube-based  antennas  to  establish  an 
effective  communication  link.  A  simulation  has  been  developed  that  estimates  the  signal-to-noise 
ratio  (SNR)  between  a  carbon  nanotube  transmitter  and  its  receiver  separated  across  a  free-space 
channel.  The  carbon  nanotube  antennas  are  modeled  using  a  thin-wire  approximation  to  the  electric 
field  integral  equation  and  solved  using  the  Method  of  Moments.  In  addition,  the  surface  impedance 
on  each  wire  is  derived  from  a  jr-spatial  tight  binding  model,  which  approximates  the  quantum 
conductive  behavior  of  the  carbon  nanotube.  To  estimate  the  Total  Channel  Loss,  the  Linville 
method  is  used  to  determine  the  power  delivered  from  the  transmit  antenna  to  the  receiver.  This 
simulation  methodology  is  used  to  perform  SNR  calculations  for  a  vertical  carbon-nanotube  dipole 
system  in  two  different  frequency  bands  (i.e.,  high  GHz,  optical).  The  performance  of  this  nano¬ 
communication  system  is  analyzed  by  conducting  a  trade-off  study  to  determine  the  required  input 
power  necessary  to  establish  an  effective  SNR. 
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1.  INTRODUCTION 


Carbon  nanotubes  show  excellent  promise  for  nanoelectronic  applications.  Because  of  small  size  and 
distinct  electrical  properties,  carbon  nanotubes  can  be  used  to  manufacture  active  device  structures, 
such  as  single-electron  transistors  [1],  or  used  as  nano-sized,  passive  structures  like  vias  and 
interconnects  [2],  [3]  and  [4],  In  this  report,  the  more  recent  application  of  exploiting  the  quasi- 
ballistic  transport  behavior  of  carbon  nanotube  for  antenna  applications  is  studied.  As  reported  in 
References  [5]  and  [6],  when  excited  electromagnetically,  a  carbon  nanotube  radiates  as  a  dipole 
antenna  according  to  its  electrical  length.  Since  the  sizes  of  carbon  nanotubes  are  easily 
manufactured  into  micron  and  nano-sized  lengths,  the  radiation  from  these  dipoles  can  be  tuned  to 
desired  frequencies  in  the  high  gigahertz  (GHz)  and  optical  frequency  bands  [8] -[10],  Research  is 
continuing  as  to  the  suitability  of  carbon  nanotubes  as  nano-antennas. 

There  have  been  several  approaches  to  studying  the  electromagnetic  properties  of  carbon  nanotube 
antennas.  Reference  [5]  developed  a  quantum  conductivity  model  for  armchair  and  zigzag  carbon 
nanotubes  and  used  the  model  to  study  the  electromagnetic  dynamics  and  surface  wave  propagation 
along  the  carbon  nanotube.  Reference  [6]  developed  a  transmission  line  model  of  a  carbon  nanotube 
antenna  and  calculated  appropriate  antenna  parameters.  In  Reference  [7],  a  finite  element  method 
was  developed  to  analyze  arrays  of  carbon  nanotube  antennas  at  terahertz  frequencies.  In  References 
[8]  through  [10],  Hanson  modeled  a  carbon  nanotube  transmitting  antenna  as  a  thin-wire  antenna.  Its 
electromagnetic  properties  were  then  determined  by  solving  a  Hallen’s  integral  equation  using  the 
Method  of  Moments.  This  approach  was  used  to  study  the  fundamental  transmitting  properties  of 
carbon  nanotubes  in  the  high  GHz  [8],  and  at  optical  frequencies  [9],  These  approaches  have 
principally  focused  on  the  transmitting  properties  of  a  carbon  nanotube  antenna.  Currently,  there  are 
no  reported  studies  that  address  the  question  of  whether  the  small  signal  strengths  that  are  being 
transmitted  by  the  carbon  nanotubes  are  sufficient  to  be  received  and  effectively  establish  a 
communication  link. 

In  order  for  a  signal  to  be  detected  and  received,  there  is  a  minimum  signal-to-noise  (SNR  ratio  that 
must  be  present  in  order  for  the  transmit  antenna  to  communicate  effectively  to  a  receiving  antenna. 
One  of  the  important  findings  of  Hanson’s  work  was  that  the  efficiencies  for  carbon  nanotube  dipole 
antennas  were  reported  to  be  very  low  (i.e.,  on  the  order  of  10'4to  10'6)  [8]— [10].  This  means  that  less 
than  .0001%  of  the  power  from  the  transmitter  is  actually  being  radiated.  In  contrast,  a  perfect 
conducting,  half-wave  dipole  will  typically  have  gains  of  approximately  2.15  dB  and  very  high 
efficiencies.  When  compared  to  the  perfect  conducting,  half-wave  dipole  antenna,  carbon  nanotube- 
based  antennas  are  clearly  poor  transmitters.  Carbon  nanotube-based  antennas  require  greater  input 
power  to  overcome  the  inherent  high  resistance  and  ohmic  losses  of  the  carbon  nanotube.  However, 
large  electric  field  strengths  may  not  be  required  to  sense,  detect,  or  communicate  at  the  nano-level. 
This  issue  is  examined  by  quantifying  the  energy  coupling  between  two  carbon  nanotube  antennas. 

The  goal  of  this  report  is  to  calculate  the  SNR  between  a  transmitter  and  receiver  antenna  system 
composed  of  two  carbon  nanotubes.  For  the  purpose  of  this  analysis,  a  communication  link  equation 
[10]  is  used  to  evaluate  the  SNR.  The  link  equation  is  expressed  as 


P  -  T 

T  T, mismatch 


Dt  +  Tjj  +  Dr  +  —  Lp  —  L 


'R, mismatch 


Nr  >  SNRrequired 


(1) 


All  values  in  Equation  (1)  are  in  dB.  Pj  is  the  signal  power  from  the  transmitter.  Lr  miSmatch  is  the 
signal  power  reflected  from  the  transmit  antenna  because  of  the  impedance  mismatch  between  the 
transmitter  and  the  transmit  antenna. 
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The  signal  power  delivered  to  the  transmit  antenna  is  \Pp  -  Lt,  mismatch ]•  Dt  and  //-/■  are  the  directivity  of 
the  transmit  antenna  and  the  transmit  antenna  efficiency,  respectively.  Similarly,  DR  and  rjR  are  the 
directivity  of  the  receive  antenna  and  the  receive  antenna  efficiency,  respectively.  LPROP  is  the 
propagation  loss  between  transmit  and  receive  antennas.  LR  miSmatch  is  the  signal  power  reflected  from 
the  receiver  load  because  of  the  impedance  mismatch  between  the  receiver  and  the  receive  antenna. 
The  Total  Channel  Loss  between  signal  power  from  the  transmitter  and  the  signal  power  delivered  to 
the  receiver  will  be  defined  as 

Total  Channel  Loss  =  Dr  +  Vt  +  DR  +  rjR  -  LPROP  -  LR  mismatch  (2) 

Finally,  Nr  is  the  receiver  noise  power.  In  order  for  a  communication  link  to  be  established,  the 
signal  power  at  the  receiver  must  be  greater  than  the  receiver  noise  power.  The  signal  power  at 
the  receiver  minus  the  receiver  noise  power  must  be  greater  than  the  required  signal-to-noise,  or 
SNRrequired.  The  SNRrequired  is  determined  by  the  properties  of  the  communication  channel  and  the 
type  of  modulation  and  demodulation  being  used. 

A  simulation  is  used  to  determine  many  of  the  link  parameters  in  Equation  (1).  The  basis  for  this 
simulation  is  the  modeling  of  the  electromagnetic  field  distribution  using  the  thin-wire  electric  field 
integral  equation  [14],  The  Method  of  Moments  [13]  is  used  to  solve  this  integral  equation.  Since 
carbon  nanotubes  are  nano-scale  devices,  the  quantum  conductive  behavior  of  the  carbon  nanotube 
must  be  incorporated  into  the  model.  An  approach,  introduced  by  Hanson,  is  used  where  the  surface 
impedance  on  the  wire  is  modeled  using  Slepyan’s  7i-spatial  tight  binding  approximation  [8].  This 
approximation  provides  an  appropriate  model  for  the  quantum  conductive  behavior  of  the  carbon 
nanotube.  Once  the  field  distribution  has  been  calculated,  the  antenna  parameters  and  signal  strengths 
can  then  be  determined  and  analyzed  for  the  communication  link.  The  electromagnetic  simulation 
and  communication  link  analysis  were  all  performed  using  a  modification  of  the  computational 
electromagnetic  code  MININEC  [18], 

The  link  equation,  Equation  (1),  will  be  evaluated  for  a  vertically  directed  carbon  nano  tube  dipole 
transmitter  and  receiver  antenna.  The  two  dipole  antennas  are  separated  across  a  free-space  channel. 
The  evaluation  is  performed  at  two  distinct  operating  ranges,  in  the  high  GHz  and  optical  frequency 
bands.  The  carbon  nanotubes  are  assumed  to  be  both  single-walled,  armchair  carbon  nanotubes. 

Based  on  this  analysis,  a  trade-off  study  is  performed  to  determine  what  input  power  is  required  to 
establish  an  effective  SNR  for  a  given  antenna  separation.  Finally,  a  study  is  performed  to  improve 
the  SNR  of  the  communication  link  by  increasing  the  directivity  of  the  carbon  nanotube  antennas.  A 
YAGI  antenna  is  considered.  This  study  is  used  is  to  estimate  the  efficacy  of  a  carbon  nanotube- 
based  antenna  communication  system. 
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2.  ELECTROMAGNETIC  MODELING  OF  CARBON  NANOTUBE 

ANTENNAS 


2.1  ELECTRIC  FIELD  INTEGRAL  EQUATION 

Figure  1  shows  a  nanotube  transmit  dipole  antenna  and  a  nanotube  receive  dipole  antenna.  The 
transmit  carbon  nanotube  is  excited  using  a  voltage  source  E'(r,  co)  =  Vs  at  the  operating  frequency  of 
co  =  27 if  This  excitation  induces  currents  on  the  nanotube  transmit  dipole  antenna.  This  current 
radiates  outward  to  generate  a  scattered  electric  field  Es(r,  co),  which  excites  and  induces  a  current  on 
the  nanotube  receive  dipole  antenna.  The  unknown  current  distribution  I(r)  induced  on  both  the 
transmitter  and  receive  carbon  nanotubes  are  calculated.  This  electromagnetically  coupled  system  is 
modeled  using  the  thin-wire  electric  field  integral  equation.  This  integral  equation  is  solved  using  the 
Method  of  Moments  [13]  for  the  current  distribution  induced  on  the  carbon  nanotubes.  Once  the 
currents  have  been  determined,  electromagnetic  properties  and  antenna  parameters  can  then  be 
computed.  The  results  can  then  be  used  to  determine  the  performance  of  a  nanotube  dipole -based 
communication  system. 


Figure  1.  Communication  system  using  carbon  nanotube  dipole  antennas. 
The  electric  field  integral  equation  for  the  coupled  antenna  system  is  given  by 

I(r) 


[E'(r,®)+E‘(r, 


(3) 


The  scattered  electric  field  Ev(r,  co)  is  expressed  in  terms  of  the  mixed  magnetic  vector  and  electric 
scalar  potentials  as 


E*  (r,  co)  =  -  jcoA  (r,w)-V$(r,ffl) 


(4) 
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In  free  space,  the  vector  and  scalar  potentials  are  defined  as 


I(r') 


-  jkR 


R 


-dS' 


(5) 


and 


<D(r  ,co)  =  ~ 


1 


jAjtcoe0 


r* 


!(■-') 
2jta[x ') 


-jkR 


R 


-dS' 


(6) 


where  R  =  |r  -  r'|  is  the  distance  between  the  source  and  observation  points.  £  =  coj/u0£0  is  the  firee- 

space  wave  number.  a( r)  is  the  radius  of  the  carbon  nanotube  and  is  assumed  constant  on  each  of  the 
two  carbon  nanotubes.  I(r)  is  the  current  vector  directed  along  the  wire  axis.  The  magnitude  of  I(r)  is 
the  total  current  on  the  wire  surface  S.  The  quantities  /u0  and  £a  are  the  permeability  and  permittivity 
of  free  space,  respectively.  By  combining  Equations  (3)  through  (6),  the  Electric  Field  Integral 
Equation  (EFIE)  for  a  carbon  nanotube-based  antenna  is  formulated  to  be 


iM 


f 


I(r') 


-jkR 


dS'  +  - 


2jra(r)  4jrJs2jta(rf)  R  j4jrco£0Js  s  \2jra(r’) 


I(r') 


\  e-M 


R 


dS’=E\m{r,co) 


(V) 


2.2  SURFACE  IMPEDANCE  BASED  ON  QUANTUM  CONDUCTIVE  BEHAVIOR 

The  surface  impedance,  Zs,  is  defined  on  the  wire,  which  takes  into  account  the  quantum  conductive 
behavior.  The  surface  impedance  Zs  in  Equation  (7)  is  derived  based  on  Slepyan’s  ^--spatial  tight 
binding  approximation  [5].  Using  this  approximation,  the  surface  impedance  Zs  per  unit  length  is 
given  by 


Z  = 


1 


2  Jta(r)oi 


CN 


(8) 


where  oqn  is  the  conductivity.  The  radius  a(r)  of  the  carbon  nanotube  is  given  by  Equation  (8), 


a 


(r)  =  — bjm 2  +  mn  +  n2 
K  ’  2. Tt 


(9) 


where  ( m ,  n)  is  the  chiral  pattern  and  b  is  the  interatomic  distance  in  grapheme,  which  is 
approximately  2.91  nm.  For  the  purpose  of  this  analysis,  armchair  carbon  nanotubes  will  be 
considered.  For  armchair  carbon  nanotubes  m  =  n.  Assuming  a  small  radii  approximation,  the 
conductance  is  given  by  the  following  formula  [5,12] 


•  2  / 

ie  co 

rr  =  — _ 

UCN  2 


Jt  ha 


-3y0b/h  ,  ,2  (Fc  -Fv)  ' 

+  2Zj*= 1  J  Ec\Rvc\  12_^_  V.A  A  J72 

-2jrh/ 


co(oj-jv) 


co(co-  jv)  -  4Ec 


Fc>v ,  Ec>y?  and  Rc>v  are  given  by  the  expressions 

Fcv=(l  +  eE",(hT)yl 


(10) 
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and 


Ecv(pz,s )  =  ±y„  ^\  +  Acquits  I  m)co?,{dp z  I  ~J3)  +  Acos1  (dp ,  /  ~J3) 


Kc(Pz’s)  =  ~oz72,  °  ,  sin(dp  /J3)sm(jrs/m) 

2Ec(pz,s) 

In  these  expressions,  e  is  the  charge  of  the  electron.  T  is  the  temperature,  v  =  1/r  is  the  relaxation 
frequency,  and  d  =  3b/2h  is  the  quasi-momentum  of  the  electron’s  conductance  in  the  z-direction. 

In  Figure  2,  the  quantum  conductivity  is  plotted  for  armchair  nanotubes  of  radii  a  =  2.71  nm  ( m  =  40) 
and  a  =  0.678  nm  ( m  =  10)  at  the  high  GHz  frequency  band.  In  Figure  3,  the  quantum  conductivity  is 
plotted  for  armchair  nanotubes  of  radii  a  =  2.71  nm  ( m  =  40)  and  a  =  0.678  nm  (m  =  10)  at  optical 
frequencies.  At  optical  frequencies,  there  are  resonant  bands  in  the  conductivity  caused  by  optical 
inter-band  transitions  found  in  the  second  term  of  Equation(lO).  The  high  GHz  frequency  band  is 
defined  as  100  to  700  GHz.  In  this  band  x  =  3.0ps.  The  optical  frequency  band  is  defined  as  400  to 
750  THz.  In  this  band,  x  =  .Olps.  For  both  of  these  frequency  bands,  the  temperature,  T,  is  295 
degrees  Kelvin  (K)  and  y0=  3.03eV. 


Figure  2.  Quantum  conductivity  for  high  GHz  of  100  to  700  GHz. 


Figure  3.  Quantum  conductivity  for  optical  frequencies  of  400  to  750  THz. 
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Figure  4  and  Figure  5  show  the  calculations  for  the  surface  impedance  resistance  (solid)  and 
reactance  (dash)  for  a  carbon  nanotube  at  optical  frequencies  (400-750  THz)  and  at  the  high  GHz 
(100-700  GHz).  The  calculations  are  for  two  different  radii,  a{r )  =  2.71nm  ( m  =  n  =  40)  in  black  and 
a{r)  =  0.678nm  (m  =  n  =  10)  in  gray.  At  the  high  GHz  of  Figure  4,  the  quantum  conductance  model 
only  depends  on  the  first  term  in  Equation  (10).  At  these  lower  frequencies,  the  surface  impedance 
does  not  depend  on  the  radius  of  the  carbon  nanotube.  Again,  at  optical  frequencies  of  Figure  5,  the 
ringing  in  the  impedances  is  caused  by  optical  inter-band  transitions. 


hT  io 


- realty 

- imagfZg) 


Z  is  independent  of  radius. 


0  100  200  300  400  500  600  700 

Frequency  (GHz) 


Figure  4.  Surface  impedance  for  high  GHz  of  100  to  700  GHz. 


Figure  5.  Surface  impedance  for  optical  frequencies  of  400  to  750  THz. 
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3.  THIN-WIRE  METHOD  OF  MOMENTS 


A  thin-wire  Method  of  Moments  is  used  to  numerically  solve  the  EFIE  in  Equation  (7)  by 
transforming  the  integral  equation  into  a  linear  system  of  equations.  The  carbon  nanotube  is  first 
discretized  into  thin-wire  segments  of  equal  radii  for  each  antenna.  The  unknown  surface  current 
across  each  wire  segment  is  expanded  in  space  using  linear  RWG  wire  basis  function  as 

l(r)=2/«A"(r)  (u) 

where  In  is  the  nth  unknown  current  coefficient  and  A„(r)  is  a  basis  function  defined  on  the  wire 
segments  [14].  Referring  to  the  geometry  in  Figure  6,  the  basis  functions  are  defined  on  a  wire 
segment  as 


A.(r)  =  yi  and  A2(r)  =  ^L  (12) 

The  quantities  D  and  §2  are  the  normalized  coordinates  of  a  wire  segment,  as  shown  in  Figure  6,  and 
are  related  by  =  1,  e  (0, 1) .  The  Jacobian  fl  is  given  by  J  =  |l '  |  for  a  wire  segment.  The  unitary 
basis  vectors  are 


1  = 


dr 

~d% 


(13) 


The  tangent  vectors  for  a  wire  segment  are  defined  in  terms  of  these  vectors  as  1 2  =  -1 1  =  1 1 .  Also,  the 
divergence  of  the  wire  basis  functions  is  given  by 


V‘A„(r)  =  l 


(14) 


Figure  6.  Element  tangent  to  a  curved  wire  segment. 

Inserting  Equation  (11)  into  Equations  (5)  and  (6)  yields  expressions  that  are  written  in  terms  of  sums 
of  partial  potentials  as 


A(r)  =  A„  (r')if  (r,r  ')do' 


(15) 


and 
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(16) 


<l>(r)  = 


-1 

jCO£ 


l 


I„fV'-A„(r')K(r,r')  da' 


where  the  Green’s  function  kernel  is  given  by 


e~jkR 

AjiR 


d(j)' 


(17) 


The  substitution  dS'  =  a(r')d<p'  do  is  made,  where  def  is  a  differential  arc  length  along  the  wire  axis  a 
and  d(f>'  is  an  angle  differential  measured  about  the  wire  circumference. 

The  vector  r'  in  Equations  (15)  through  (17)  is  reinterpreted  as  a  vector  to  a  point  on  the  wire  axis. 
This  means  that  R  is  no  longer  |r  -  r  j,  but  is  the  distance  from  r  to  a  point  on  the  wire’s  surface  with 
its  circular  cross  section  centered  at  r\  Also,  note  that  geometrical  complexities  such  as  bends,  radii 
discontinuities,  and  junctions  that  can  arise  at  wire  segment  endpoints  are  ignored;  it  is  merely 
assumed  that  total  current  is  continuous  (or,  equivalently,  that  Kirchhoff  s  current  law  is  satisfied  at  a 
junction).  There  are  no  issues  with  using  Equations  (15)  through  (17)  when  the  source  and 
observation  points  are  not  close.  However,  when  the  observation  point  approaches  the  source  point, 
the  wire  kernel  becomes  singular.  In  the  past,  the  singular  part  was  extracted  in  some  way  and 
evaluated  analytically.  With  the  recent  transformation  presented  in  References  [14]  and  [15], 
singularity  extraction  is  no  longer  needed  since  the  inner  integral’s  integrand  is  smoothed  by  variable 
transformations,  whereas  the  logarithmically  singular  outer  integral  may  be  integrated  numerically 
using  the  scheme  described  in  Reference  [15]. 

Next,  a  suitable  testing  scheme  must  be  adopted.  In  this  formulation,  standard  Galerkin  testing 
functions  are  used  in  space,  evaluated  via  a  higher-order  Gaussian  rule.  Applying  Galerkin  testing  to 
the  EFIE  in  Equation  (7)  yields  the  following  linear  system  of  equations, 


2  ^*4^  (*~  )  ,/y  Up  r  pAm(r)g\B  (**  ) 

4  2jia(r)  4 2jra(r') 

,  1  rrVs-A„(r)V'-l(rV« 

j4jrco£0  2  jia  ( r ' )  R 


e-jkR 

- - dS’dS 

R 


dS'dS  -ia-  (r)'E*tan  (r)dS',  m  =  l,2,K,7V 


(18) 


Finally,  the  partial  potentials  of  Equations  (15)  through  (17)  are  substituted  into  the  Galerkin  tested 
EFIE  given  by  Equation  (18).  This  yields  the  following  linear  system  of  equations, 

^ m,n I n  =  ?  TVl  =  1,  2,K  9N  (19) 


which  can  be  solved  to  determine  the  unknown  current  coefficients  In  located  at  the  nodes  between 
two  wire  segments  making  up  each  antenna.  The  matrix  ZmjU  in  Equation  (19)  is  an  TVx  N  system  of 
equations  given  by 


Am(r)gV„(r)  n  Am(r)g\„  (r')  e-*« 
SJS  2jza(r')  2jta(r')R  R 


dS'dS 


+ 


1  rrVs- Am(r)V' -A„(r')e-^ 
4  2jtci[y')R  R 


dS'dS , 


/77,/7  =  1,2,K  ,N 


(20) 


The  right-hand  side  voltage  Vm  in  Equation  (18)  is  used  to  excite  and  is  written  as 
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(21) 


Fm-fAm(r)-El(r)dS 

s 

For  this  analysis,  the  transmit  antenna  is  excited  using  a  voltage  delta-gap  generator  source.  In  this 
case,  the  wire  segments  attached  to  a  current  node  (junction)  m  can  be  thought  of  as  separated  by  a 
gap  across  which  the  voltage  Vm  is  specified. 
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4.  FUNDAMENTAL  ANTENNA  PROPERTIES 
OF  CARBON  NANOTUBES 


4.1  IMPEDANCE 

The  impedance  of  a  nanotube  transmit  dipole  antenna  is  shown  in  Figure  7.  The  half-length  of  the 
dipole  antenna  is  10  pm.  The  radius  is  2.71  nm.  This  radius  corresponds  to  an  armchair  carbon 
nanotube  where  m  =  n  =  40.  The  resistance  varies  from  approximately  10  to  140  kilo-ohms  (i.e., 
thousands  of  ohms).  The  reactance  varies  from  -80  to  50  kilo-ohms. 


Figure  7.  Impedance,  high  GHz. 

4.2  FAR-FIELD  RADIATION  PATTERN 

Once  the  induced  currents  on  the  carbon  nanotubes  antennas  have  been  computed,  the 
electromagnetic  properties  can  be  analyzed.  The  far-field  electric  field  is  approximated  as 

E(R)  =  "  /(5)]  eJka'dS  (22) 

where  k  =  (cos0  cosip)  x  +  (cosO  simp)  y  +  (simp)  z  is  the  wave  vector,  tj  =  a lju0/e0  is  the  free-space 
impedance.  The  vector  R  is  the  range  vector  from  the  current  on  the  carbon  nanotube  wire  to  the  far- 
field  observation  point.  The  integral  is  evaluated  in  closed  form  over  each  straight  wire  segment.  The 
total  far-field  radiation  pattern  is  then  the  sum  of  the  contributions  from  all  of  the  wire  segments. 
Based  on  the  far-field  electric  field  calculation,  the  radiation  pattern  of  the  antenna  is  given  by 

P(0,<p)  =  ±R2RQ[ExH]  =  ^Eg:*  (23) 

In  Figure  8,  the  current  distribution  on  a  carbon  nanotube  dipole  antenna  is  plotted  at  200,  400,  and 
600  GHz.  Again,  the  half-length  of  the  dipole  antenna  is  10  pm  and  has  a  radius  of  2.71  nm  (m  =  n  = 
40).  Figure  9  shows  the  corresponding  far-field  radiation  pattern  for  the  dipole  at  these  frequencies. 
All  of  the  far-field  patterns  of  Figure  9  show  the  classic  sin  (0)  pattern  distribution  of  the  small 
dipole. 
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Figure  8.  Current  distribution,  10-gm  half-length  dipole  antenna  of  radius  of  2.71  nm. 


±180° 

-165°  -0.4x1 02 165° 


Figure  9.  Far-field  pattern,  10-pm  half-length  dipole  antenna  of  radius  of  2.71  nm. 
The  three-dimensional  radiation  pattern  of  the  dipole  at  200  GHz  is  shown  in  Figure  10. 
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Figure  10.  3-D  far-field  pattern,  10-gm  half-length  dipole  antenna  of  radius  of  2.71  nm. 

4.3  ELECTRIC  NEAR  FIELD 

To  provide  better  understanding  of  the  magnitude  of  field  strengths,  the  electric  near  fields  in  the 
vicinity  of  the  antenna  are  simulated.  The  electric  near  fields  can  be  determined  from  the  current 
distribution  by  using  the  virtual  dipole  method  [17].  To  determine  the  electric  field  at  a  given  point  in 
the  near  field,  a  small  virtual  thin-wire  dipole  is  placed  at  the  point  with  its  axis  parallel  to  the 
appropriate  vector  component.  The  open-circuit  voltage  at  the  near-field  point  can  be  calculated  from 
knowledge  of  the  calculated  current  distribution  and  the  mutual  impedance  between  the  carbon 
nanotube  antenna  and  the  virtual  dipole.  In  other  words, 


(24) 


The  virtual  dipole  is  open-circuited,  and  Vd  is  the  open-circuit  voltage.  I,  designates  the  computed 
currents  on  the  wire  structure.  Zdi  designates  the  mutual  impedances  between  the  wire  structure  and 
the  virtual  dipole.  The  mutual  impedances  are  calculated  from  the  Method  of  Moments.  The  electric 
field  strength  along  the  direction  of  the  virtual  dipole  of  length  Ld  is 

Ed=(Vd/Ld) 

The  equation  is  evaluated  once  for  each  electric  field  component  in  the  x,  y,  and  z  directions  at  the 
near-field  point  of  interest.  The  virtual  dipole  length  is  .01  of  the  shortest  segment  length  in  the 
antenna  segmentation. 

The  three  electric  field  vector  components,  Ex ,  Ey,  and  EZi  are  calculated  as  real  and  imaginary  terms, 
from  which  the  magnitude  and  phase  are  determined.  The  average  value  of  the  electric  field  is 
defined  by  Eave  =  \E\.  The  maximum  or  peak  electric  field  is  [17] 


u  peak 


±(E2+Ey2+E2)  +  ±(A2+B2) 


'2x1/2 


1/2 


(25) 


where 

A  =  E2x  cos(2 0x)  +  E2  cos(2 6y)  +  E2  cos(26z) 
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and 


B  =  El  sin(2/9x)  +  Ey  sin(2^)  +  Ez  sin(2 Oz) 

and  where  6X ,  9y ,  and  0Z  are  the  phase  angles  for  the  corresponding  field  components. 

In  Figure  1 1  and  Figure  12,  the  electric  near  field  (x  and  z-components)  of  the  10-pm  half-length 
nanotube  dipole  antenna  and  radius  of  2.71  nm  are  shown.  The  dipole  has  been  excited  by  1  watt  of 
power  at  200  GHz. 
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Figure  11.  Electric  near  field  Exin  dB,  10-gm  half-length  dipole  antenna  of  radius  of  2.71  nm. 
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Figure  12.  Electric  near  field  £zin  dB,  10-pm  half-length  dipole  antenna  of  radius  of  2.71  nm. 

4.4  EFFICIENCY 

There  are  two  methods  for  calculating  the  efficiency  of  a  transmitting  antenna  structure.  In  the  first 
method,  the  far-field  power  gain  is  calculated.  The  far-field  power  gain,  in  a  given  direction  in 
spherical  coordinates,  is 
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(26) 


Gdb(o,(p)=mogm 


4jrP(ft,(p) 

4 


where  G^{6,  <?>)  is  the  power  radiated  per  unit  steradian  in  the  angular  directions  6  and  cp.  Pin  is  the 
total  input  power  to  the  antenna  calculated  from  the  applied  voltages  Vs  and  the  corresponding  feed 
point  currents  I  s  as 


/>„-Re(Fs/;)/2  (27) 

Directive  gain  of  a  transmit  antenna  is  the  power  gain  divided  by  the  efficiency  of  the  antenna,  or 

4/s  =  4/s  -  JldB  (28) 

where  ridBis  the  efficiency  in  dB.  The  directive  gain  of  a  small  antenna  is  1.77  dB.  The  difference 
between  this  directive  gain  and  the  calculated  power  gain  would  be  the  efficiency. 

In  the  second  method,  the  ohmic  power  losses  on  the  antenna  structure  can  be  calculated.  The  ohmic 
power  loss  is  calculated  from  the  currents  of  Equation  (7)  multiplied  by  the  real  part  of  the  surface 
impedance,  using 


4TO-c=l/2|42Re(4)  (29) 

The  input  impedance  of  a  carbon  nanotube  antenna  is  given  by  Zin  =  Vs  /  Is.  The  efficiency  is  then  the 
radiated  power,  or  the  input  power  minus  the  ohmic  power  loss  divided  by  the  input  power,  or 

p  ~p 

V  -  —  (30) 

*in 

The  two  methods  of  analysis  yielded  similar  efficiency  values. 

4.5  IMPEDANCE  MISMATCH  LOSS 

Another  loss  is  associated  with  the  reflected  power  back  to  the  transmitter  because  of  impedance 
matching  between  the  antenna  impedance  and  the  internal  impedance  of  the  transmitter.  This 
mismatch  loss  reduces  the  power  delivered  to  the  transmit  antenna.  When  a  source  is  terminated  with 
an  antenna  with  an  impedance  of  ZAT ,  that  is  not  equal  to  the  internal  impedance  of  the  transmitter  or 
source,  ZT ,  not  all  of  the  power  is  delivered  to  the  antenna.  This  reflected  power  results  in  a  mismatch 
loss  of 


^ Tmismatch  201og 


7  -7 

^AT 

Zat  +  ZT 


(31) 


The  input  impedance  is  the  source  voltage  divided  by  the  current  flowing  at  the  source.  For  the 
purpose  of  this  analysis,  ZT  =  12,900  ohms.  This  impedance  value  is  based  on  the  nanoscale,  dc 
electron  transport  [6] .  A  similar  mismatch  can  be  defined  for  the  mismatch  between  the  receive 
antenna  impedance  and  the  represented  load  of  the  receiver.  This  receive  mismatch  loss  is 
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^ Rmismatch  201Og 


(32) 


^AR 

Z AR  +  Zr 


4.6  RECEIVER  NOISE  POWER 

Finally,  for  the  purpose  of  this  analysis,  it  is  assumed  that  the  receiver  noise  power  is  generated  in  the 
receiver  itself.  The  contribution  of  external  noise  is  minimal.  The  receiver  noise  power,  Nr,  is  given 
by 


Nr  -kTRB  (33) 

where  k  is  Boltzmann’s  constant.  Tr  is  the  receive  system  temperature  and  B  is  the  receiver 
bandwidth.  The  receive  noise  temperature  is  usually  determined  from  the  Noise  Figure  in  dB,  NFdB, 
of  the  receiver  and  the  ambient  temperature,  290  degrees  Kelvin.  The  receive  system  temperature  can 
be  determined  from 


Tr  =  290 


NFm\ 

10 


-1 


(34) 


4.7  TOTAL  CHANNEL  LOSS 


In  Equation  (2),  the  Total  Channel  Loss  was  defined  as  [Dj+  >]t+  Dr  +  >]r  -  Lprop  -  Lr,  mismatch]- 
Rather  than  calculating  each  parameter  in  the  Total  Channel  Loss,  it  is  possible  to  calculate  the  Total 
Channel  Loss  directly  using  the  Linville  method  [11]. 

The  first  step  is  to  calculate  the  admittance  parameters  for  the  transmitter  and  receiver  antenna 
system.  Any  two  antennas  can  be  treated  as  a  two-port  network,  where  terminals  1-1'  would  be  the 
feed  point  of  the  first  antenna  and  the  terminals  2-2'  would  be  the  feed  point  of  the  second  antenna. 
The  two-port  network  equations  can  be  written  as 


m 

[7n 

Kll 

722. 

(35) 


These  admittance  parameters  are  defined  by  the  following  relations 


Y  =- LI 
11 


Vj=0 


y12=L 
l!  V, 


K=0 


v  - 1 2 1 

y} 


F,  =0 


7  -h, 

’  22  “  K 


F,=0 


where  I  )  and  Vi  equal  zero  implies  that  the  terminals  associated  with  these  voltages  are  short- 
circuited.  The  current  calculation  is  used  to  find  the  admittance  parameters.  The  feed  point  of  one 
antenna  is  excited,  and  the  feed  point  of  the  second  antenna  is  short-circuited.  The  calculated  currents 
are  used  to  determine  Yu,  Yu,  7 21,  and  722.  Once  the  admittance  parameters  have  been  found,  the 
Total  Channel  Loss  can  be  calculated  for  arbitrary  receiver  loads. 

To  evaluate  the  maximum  power  transfer  between  the  transmitter  and  the  power  received  at  the 
receiver  or  minimum  Total  Channel  Loss,  the  input  and  output  power  must  be  determined.  Given  an 
unknown  load  on  the  receive  antenna,  the  output  power  is  determined  by 
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(36) 


2 

p  = 

out 

Vr 

Re 

^load 

where  Vr  is  the  output  voltage  on  the  receive  load,  and  Re  \Yioad  ]  is  the  real  part  of  the  receive  load 
admittance  Fioaci.  Similarly,  the  input  power  can  be  calculated  by 

^=|^r|2Re[4]  (37) 

where  Vj  is  the  input  voltage  on  the  transmit  antenna.  Re  [  is  the  real  part  of  the  input  admittance 
of  the  transmit  antenna  port.  Taking  the  ratio  of  Equations  (36)  and  (37),  a  channel  gain  can  be 
defined  as 


3 Channel 


R 


Vr 

Re  [4] 

vT 

This  ratio  is  dependent  on  the  input  admittance  Yin  given  by 


Y  Y 

V  _V  1 RT1TR 

1in  ~  1TT  y  ,y 

1  load  'riRR 


(38) 


(39) 


The  two-port  admittance  parameters  can  then  be  used  also  to  find  the  output  to  input  voltage  ratio,  Vr 

/  Ft,  as 


VR 

Y  Y 

1 RT1TR 

VT 

Y  +Y 

1  load  ^  1 RR 

(40) 


Substituting  Equations  (39)  and  (40)  into  Equation  (38),  the  channel  gain  becomes 


J  Channel 


Ref  M 

V  y 

1 RT1TR 

Re  [4] 

^load  +  YrR 

(41) 


In  order  to  minimize  the  Total  Channel  Loss,  the  load  admittance,  7/oaj,  must  be  calculated  that 
maximizes  the  channel  gain.  The  maximum  channel  gain  is  given  by 


where 


L  = 


P 

ll 

r  /  o\1/2l 

l-(l-L2) 

If  y  1 

\1RT1TR\ 

2Re  Ytt  Re  YRR  -Re\YRTYTR 


(42) 


(43) 
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5.  NUMERICAL  RESULTS 


The  performance  of  a  communication  system  using  identical  carbon  nanotube  antennas  is  analyzed 
using  the  equations  from  the  above  sections.  Results  are  presented  using  vertical  dipole  antennas  and 
YAGI  antennas.  The  antenna  systems  are  operated  in  both  the  high  GHz  and  at  optical  frequencies. 

5.1  DIPOLES  IN  THE  HIGH  GHZ 

Results  are  presented  for  a  communication  system  using  parallel  transmit  and  receive  carbon 
nanotube  dipole  antennas,  both  with  half-lengths  of  10  pm.  Both  transmit  and  receive  antennas  are 
assumed  to  be  single-walled,  armchair  carbon  nanotubes  with  a  radius  of  10  nm  (m  =  n  =  50).  The 
transmitter  and  receiver  antennas  are  separated  in  free  space  across  a  channel  of  length  R  that  varies 
from  10  pm  to  10,000  pm  in  distance.  The  transmit  antenna  is  excited  at  the  center  of  the  antenna. 
The  receive  dipole  is  loaded  at  the  center  of  the  dipole.  The  frequency  range  is  in  the  high  GHz 
regime  from  100  GHz  to  1000  GHz. 

In  Figure  13,  the  normalized  impedance  (Z„,  /  R„)  for  the  carbon  nanotube  transmitter  antenna  is 
shown.  The  left  and  right  axis  of  the  plot  is  for  the  real  part  and  imaginary  part  of  the  normalized 
impedance,  respectively.  The  resonance  of  the  dipole  is  determined  when  the  imaginary  part  of  the 
impedance  crosses  zero.  The  dipole’s  first  resonance  and  anti-resonance  occurs  at  290  GHz  and  600 
GHz,  respectively.  These  values  are  comparable  to  those  results  presented  by  Hanson’s  solution  to 
Hallen’s  integral  for  the  same  antenna  [8],  Note,  for  many  of  the  calculations  below,  the  impedance 
is  normalized  using  Ra,  12,900  ohms.  Again,  this  normalized  impedance  value  is  based  on  the 
nanoscale,  dc  electron  transport  [6]. 


Frequency  (GHz) 


Figure  13.  Normalized  impedance  of  carbon  nanotube  transmitter  dipole,  high  GHz. 

Using  Equation  (29),  the  ohmic  losses  on  the  transmit  antenna  can  be  determined  from  the  Zs  and  the 
computed  currents.  In  addition,  the  mismatch  loss  between  the  transmitter  and  the  transmit  antenna 
can  also  be  determined  from  Equation  (31).  ZT  is  assumed  to  be  12,900  ohms.  Both  of  these  results 
are  plotted  in  Figure  14.  The  ohmic  losses  are  considerably  higher  than  the  mismatch  loss.  The 
carbon  nanotube  antenna  is  fairly  well  matched  to  the  internal  impedance  of  a  nano-transmitter  with 
an  internal  impedance  of  12,900  ohms.  On  the  other  hand,  the  high  surface  impedance  of  a  carbon 
nanotube  antenna  is  a  significant  source  of  power  loss. 
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Figure  14.  Ohmic  and  mismatch  loss  of  carbon  nanotube  transmit  dipole,  high  GHz. 

Since  the  ohmic  losses  on  the  carbon  nanotube  are  so  high,  the  radiation  efficiency  of  the  antenna  is 
expected  to  be  very  poor.  In  Figure  15,  the  efficiency  for  the  transmit  dipole  antenna  is  plotted  as  a 
function  of  frequency.  The  dipole’s  efficiency  has  a  minimum  of  about  -93  dB  at  its  first  anti¬ 
resonance,  which  occurs  around  600  GHz.  Its  maximum  efficiency  is  at  the  first  resonance,  290  GHz, 
at  a  calculated  value  of  -50  dB.  This  means  that  only  .000025  %  of  the  input  power  is  actually  being 
radiated.  Most  of  the  input  power  is  absorbed  by  ohmic  loss  on  the  carbon  nanotube. 


Frequency  (GHz) 


Figure  15.  Efficiency  of  carbon  nanotube  transmit  dipole,  high  GHz. 

The  Total  Channel  Loss  of  Equation  (2)  is  shown  in  Figure  16.  The  transmitter  and  receiver  are 
separated  by  distances  varying  from  10  pm,  which  is  in  the  near  field  of  the  antenna,  to  1,000  pm.  At 
the  first  resonance  at  290  GHz,  the  Total  Channel  Loss  is  nominally  varying  from  between  -50  to 
-150  dB  as  the  distance  between  the  antennas  increases.  Because  the  efficiency  of  the  transmit 
antenna  is  so  poor  at  its  first  anti-resonance  of  600  GHz,  there  is  very  little  power  that  is  transmitted 
across  the  channel  and  the  Total  Channel  Loss  varies  from  between  -50  to  -225  dB.  This  Total 
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Channel  Loss  is  for  the  optimal  receiver  load.  However,  there  was  minimal  difference  between  the 
Total  Channel  Loss  for  the  optimal  receiver  load  and  a  receiver  load  of  12,900  ohms. 


Figure  16.  Total  Channel  Loss  between  transmit  and  receive  carbon  nanotube 

dipole  antennas,  high  GHz. 

Again,  considering  Equation  (1),  a  communication  link  analysis  is  performed.  For  the  purpose  of  this 
analysis,  a  Noise  Figure  of  7  dB  is  assumed.  The  required  signal-to-noise  ratio  is  assumed  to  be  5  dB. 
These  are  nominal  values  for  the  purpose  of  this  study.  Both  values  are  optimistic  for  a  nano- 
technology-based  communication  link.  Using  Equation  (1),  the  required  transmit  or  input  power  in 
watts  can  be  determined  for  the  specified  parameters.  The  results  are  displayed  in  Figure  17.  At  600 
GHz  and  a  distance  of  1000  pm  between  the  transmitter  and  receiver,  approximately  105  watts  of 
input  power  is  required  for  effective  communications.  It  is  apparent  that  the  ohmic  losses  severely 
limit  the  ability  to  effectively  communication  between  two  nanotube  dipoles. 


Figure  17.  Required  input  power  between  transmit  and  receive  carbon  nanotube 

dipole  antennas,  high  GHz. 


19 


5.2  VERTICAL  DIPOLES  AT  OPTICAL  FREQUENCIES 


Results  are  presented  for  the  vertically  directed  dipole  antenna  systems  operating  at  optical 
frequencies  (i.e.,  100  THz  to  800  THz).  Both  carbon  nanotube  dipole  antennas  have  half-lengths  of 
150  nm  with  a  radius  of  2.72  nm  ( m  =  n  =  40).  The  transmitter  and  receiver  antennas  are  separated  in 
free  space  across  a  channel  of  length  R  that  varies  between  100  nm  to  10,000  nm  in  length.  Again, 
the  transmit  antenna  is  excited  at  the  center  of  the  antenna.  The  receive  dipole  is  loaded  at  the  center 
of  the  dipole.  Figure  18  shows  the  normalized  impedance  for  the  transmitting  dipole  antenna.  At 
these  frequencies,  the  dipole  antenna  is  not  resonating.  The  ringing  in  the  impedance  is  caused  by  the 
conductivity  of  the  carbon  nanotube  resonating  in  the  optical  bands. 

Figure  19  shows  the  efficiency  of  the  transmit  dipole  antenna  at  optical  frequencies.  The  carbon 
nanotube  antenna  still  remains  very  inefficient  at  these  frequencies. 


Figure  18.  Normalized  impedance  of  carbon  nanotube  transmit  dipole,  optical. 


Figure  19.  Efficiency  of  carbon  nanotube  transmit  dipole,  optical. 
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The  Total  Channel  Loss  is  shown  in  Figure  20.  The  transmitter  and  receiver  are  separated  by 
distances  L  =  100  nm,  1000  nm,  and  10,000  nm.  At  120  THz,  the  Total  Channel  Loss  is  nominally 
varying  from  between  -80  to  -150  dB  as  the  distance  between  the  antennas  increases. 


Figure  20.  Total  Channel  Loss  between  transmit  and  receive  carbon  nanotube 

dipole  antennas,  optical. 

5.3  ANTENNAS  OF  HIGHER  DIRECTIVITY 

A  final  study  is  performed  to  improve  the  SNR  of  the  communication  link  by  using  antennas  that 
would  have  higher  directivity  than  the  dipole  antenna.  As  reported  in  [8],  antennas  that  are 
considered  to  have  higher  directivity  were  expected  to  increase  the  magnitude  of  the  electric  field 
strength.  YAGI  antennas  have  higher  directivity  than  the  dipole  antenna.  A  comparison  is  performed 
between  the  nanotube  dipole  and  a  three-element  nanotube  YAGI  antenna.  The  specific  three- 
element  YAGI  antenna  would  have  a  directive  gain  of  approximately  9.7  dB.  This  would  occur  if  the 
YAGI  antenna  is  metal  and  operating  at  its  first  resonance. 

Figure  21  shows  the  plotted  efficiencies  for  the  two  different  antennas.  Surprisingly,  there  is  very 
little  variation  between  the  antennas.  For  the  two  antennas,  the  maximum  efficiency  is  still  around 
290  GHz,  at  a  calculated  value  of  -50  dB. 
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Figure  21.  Efficiencies  of  dipole  (bold)  and  YAGI  (+)antennas,  high  GHz. 

Finally,  a  communication  link  analysis  is  performed  for  the  two  antennas.  The  Total  Channel  Loss  is 
for  the  optimal  receiver  load.  Again,  a  noise  figure  of  7  dB  is  assumed.  The  required  SNR  is  assumed 
to  be  5  dB.  The  transmitter  and  receiver  antennas  are  separated  by  distances  of  100  pm,  1000  pm, 
and  10,000  pm.  The  required  input  power  in  watts  for  the  specified  parameters  is  shown  in  Figure  22. 
There  appears  to  be  no  variation  among  the  two  antennas.  It  is  apparent  that  the  ohmic  losses 
severely  limit  the  transmitted  signal.  Antennas  of  higher  directivity  have  limited  utility  in  improving 
the  communication  link  performance. 


Figure  22.  Required  input  power  for  dipole  (bold),  and  YAGI  (+)  antennas,  high  GHz. 


22 


6.  CONCLUSIONS 


A  simulation-based  study  has  been  presented  for  the  communication  link  between  carbon  nanotube- 
based  antennas.  The  approach  is  based  on  calculating  a  link  equation  given  the  transmitter  and 
receiver  antenna  properties  and  losses  generated  in  the  system.  Using  the  link  equation,  various 
feasibility  studies  can  be  performed,  including  estimating  signal-to-noise  ratios  and  determining 
minimum  discernable  signals.  However,  the  focus  of  this  study  was  on  estimating  the  required  input 
power  necessary  to  establish  a  communication  link  given  a  SNR  greater  than  zero.  Because  of  the 
high  ohmic  losses  generated  on  the  antennas,  considerable  and  highly  improbable  input  power  is 
necessary  to  establish  the  link.  In  order  to  maintain  a  SNR  of  5  dB,  input  power  on  the  order  of  10  to 
le5  watts  was  required.  While  it  is  believed  that  the  performance  of  the  communication  system  could 
be  further  improved,  it  is  highly  unlikely  that  any  optimized  configuration  could  overcome  the  high 
ohmic  losses. 
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